Ground State in a Half-Doped Manganite Distinguished by Neutron Spectroscopy 
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We have measured the spin wave spectrum of the half-doped bilayer manganite Pr(Ca,Sr)2Mn207 
in its spin, charge, and orbital ordered phase. The measurements, which extend throughout the 
Brillouin zone and cover the entire one-magnon spectrum, are compared critically with spin-wave 
calculations for different models of the electronic ground state. The data are described very well by 
the Goodenough model, which has weakly interacting ferromagnetic zig-zag chains in the CE-type 
arrangement. A model that allows ferromagnetic dimers to form within the zigzags is inconsistent 
with the data. The analysis conclusively rules out the strongly bound dimer (Zener polaron) model. 

PACS numbers: 75.30.Ds, 71.10.-w, 75.25.Dk, 75.47.Lx 



Perovskite manganese oxides and their layered ana- 
logues continue to be widely studied not only for 
the colossal magnetoresistance (CMR) they can show 
(changes in resistivity of up to several orders of magni- 
tude in magnetic fields of a few Tesla) but also because of 
the rich physics arising from the coupling of spin, charge, 
orbital, and lattice degrees of freedom pj. The most pro- 
nounced CMR phenomena typically arise when the man- 
ganite is delicately poised between a ferromagnetic metal 
and an antiferromagnetic (AFM) charge-orbital (CO) or- 
dered insulator. One of the most ubiquitous AFM CO 
phases is found in manganites near half doping Q , yet de- 
spite extensive study motivated by the link with CMR, 
the electronic state and the mechanisms that stabilize 
this AFM CO phase remain unclear. 

The first attempt to describe the AFM CO phase was 
by Goodenough, whose celebrated model is depicted in 
Fig. [Ija). It is based on the Goodenough-Kanamori- 
Anderson rules for superexchange 0, 0] and assumes a 
checkerboard pattern of Mn'^^ and Mn''^"'" charge order 
, and an associated herringbone pattern of occupied Cg 
■ orbitals on Mn'^"'" sites. This results in the CE-type 
^ , magnetic order, comprising ferromagnetic (FM) zig-zag 
chains of spins with antiferromagnetic (AFM) alignment 
to each neighboring chain. Neutron diffraction [5|-|7| and 
resonant magnetic x-ray scattering studies appear 

to support this model, but with markedly reduced (to 
< 25%) charge disproportionation 0, 0, El-Eil ■ ^ dis- 
tinctly different model, the Zener polaron (ZP) model 
[iHllali postulates that pairs of nearest- neighbor Mn ions 
effectively share an electron through the Zener double- 
exchange mechanism. These polarons form the bond- 
centred charge ordering pattern shown in Fig. [Ijb) , the 
spins within a polaron being strongly coupled ferromag- 
netically |16l4l9| . While some structural studies come 

others 



out strongly against the ZP model [12|, [1 



support mixed valence or the ZP model 



plex ordering periodicities, and are very difficult to dis- 
tinguish in diffraction experiments. An alternative ap- 
proach is to study the spectrum of magnetic excitations, 
which is determined by the strength of the exchange in- 
teractions and will therefore reflect patterns of orbital 
ordering and dimerisation. In particular, the CE-type 
and ZP magnetic structures have very different spin- 
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The Goodenough and ZP models have the same corn- 



wave spectra [24| . Previous neutron scattering mea- 
surements of the spin-wave dispersion in Lao.5Sri.5Mn04 
are in good agreement with the Goodenough model [2^ . 
However, these measurements neglect the higher energy 
(> 40meV) part of the spin- wave dispersion, and cannot 
conclusively rule out the possibility of a Zener polaron 
state (sec below). The same is true of low energy spin- 
wave data on Ndo.sSro.sMnOg [2^. 

Here we report neutron scattering measurements of the 
complete magnetic spectrum of a half-doped bilayer man- 
ganite. We compare the results quantitatively with spin- 
wave calculations for the Goodenough and ZP models. 
The results provide a clear preference for the Goode- 
nough model (but with reduced charge disproportiona- 
tion) over both the ZP model and an alternative model 
with weaker dimerisation in the FM zigzags. 

The AFM CO phase considered here is observed at 
or near half doping in the Ruddlesden-Popper series 
with general formula, A„_|_iMn„03„+i, where ^ is a rare- 
earth-metal or alkaline-earth-metal ion. The n = 2 mem- 
ber of the series, considered here, has a bilayered struc- 
ture [see Fig. Hfd)]. The magnetic coupling between 
Mn204 bilayers is expected to be < 1/100 of the intra- 
bilayer interactions [28|. As well as having a quasi- two- 
dimensional (quasi-2D) magnetic dispersion, the layered 
systems have only two crystal twins rather than six as 
in the three dimensional (3D) manganites (n = 00) [l5| . 
which greatly simplifies data analysis. 

Our experiment was performed on half-doped 
Pr(Cao.9Sro.i)2Mn207 (PCSMO). PCSMO exhibits 
charge and orbital order below Tqq = 370 K and antifcr- 
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FIG. 1. (Color online). Magnetic order and exchange inter- 
actions within a Mn-0 layer for the three models considered 
here: (a) the Goodenough model with CE-type magnetic or- 
der, (b) the Zener polaron model, and (c) the dimer model. 
The in-plane magnetic unit cells are indicated by the dotted 
lines, (d) Crystal structure of PCSMO [H, with Mn ions 
in oxygen octahedra (green) and the Pr, Ca, or Sr ions (yel- 
low) between the Mn-O layers, (e) Projection of the high 
temperature Amam unit cell onto the Mn-O layer. 



romagnetic order below Tn = 153 K 2^31|. The specific 
composition of PCSMO used here can be prepared in 
the form of large crystals suitable for inelastic neutron 
scattering. 

Neutron diffraction data confirmed the presence in our 
sample of two sets of magnetic Bragg peaks below Tn 
characteristic of the AFM CO phase magnetic order, one 
with indices (0.5,0,/), and the other (0.5,0.5,/), with / 
an integer Q. Throughout this Letter we use the high 
temperature orthorhombic Amam unit cell (29l |. with 
a = 5.41 A, b = 5.46 A, and c = 19.3 A. In this setting, 
the a axis is parallel to the CE-type FM zigzags, [see 
Fig. [ije)]. Because of twinning, magnetic peaks are also 
observed at equivalent positions with the a and b axes 
interchanged. The diffraction data also showed that the 
magnetic coupling between layers within a bilayer (J_l) 
is AFM. 
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FIG. 2. (Color online). Neutron scattering intensity maps 
of the spin-wave dispersion of PCSMO in the {h, 0) direc- 
tion. The left-hand plots show energies up to 30meV, and 
the right-hand plots show the whole dispersion. Panels (a) 
and (b) are experimental data, and panels (c) and (d) are 
calculations from spin-wave theory for the CE-type order (see 
text). Regions of the spectrum demarcated by dashed lines 
were measured with different neutron incident energies and 
(in some cases) different Brillouin zones. The data have been 
averaged over a small range of wave vectors in the (0, k) 
direction. To display the full spectrum with one intensity 
scale, we have multiplied the data by the monotonic func- 
tion f{E) = Yzf^WW' Eo = kBX 10 K, which has the effect 
of suppressing the strong elastic and low energy scattering. 
Intensity units are mbsr"^ meV~^ per Mn. 



Neutron scattering spectra were recorded at a temper- 
ature of 4 K on a time-of- flight spectrometer in a configu- 
ration that made it possible to measure the in-plane mag- 
netic dispersion as a function of neutron energy trans- 
fer (see Supplemental Material [s^). Figures [IJ^a) and 
2(b) present intensity maps as a function of energy and 
wave vector that summarise the important features of 
the spin waves. Because the dispersion is quasi- 2D we 
give only the in-plane wave vector components. The data 
have been corrected for an energy-dependent background 
taken from positions where magnetic scattering is negli- 
gible. The magnetic dispersion is plotted in the direction 
parallel to the CE-type FM zigzags, but also contains 
features dispersing perpendicular to the zigzags due to 
twinning. There is a low energy band emerging from the 
magnetic Bragg peak at (0.5, 0) with a maximum energy 
of 7meV at (0,0), and an upper band which disperses 
from 7meV at (0,0) to a maximum energy of 85meV 
at (1,0). The upper band has an energy gap between 
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FIG. 3. (Color online). Constant wave vector cuts through 
the experimental data. Fits with resolution corrections are 
shown by the (blue) lines. A non-magnetic background has 
been subtracted. 



35 and 55 mcV, which extends throughout the Brillouin 
zone (not all shown). The low energy band is from the 
direction perpendicular to the FM zigzag, and all of the 
dispersion over 7meV is from the direction parallel the 
FM zigzag. The periodicity of the upper band, Ah = 2, 
immediately reveals that the dominant exchange interac- 
tion is FM between nearest-neighbor Mn ions projected 
along the {h, 0) direction. The presence of the gap at 
h = 0.5 requires a real-space doubling of the period, for 
example from dimerisation of Mn ions, second-neighbor 
interactions, or a different moment on alternate Mn ions. 
A quasi-one-dimensional (quasi-lD) magnetic exchange 
is implied by the fact that for all energies above '^12 
meV, the intensity in {h,k) planes at constant energy 
can be explained by dispersion independent of k once 
the twinning is accounted for. 

Figure [3] shows some examples of energy cuts through 
the data at different wavevectors. Figures [31Ja)-[31[c) con- 
tain data on the lower band, and Figs. [31[d)-[3][f) on the 
upper band. The cut in Fig. ^a) shows that there is 
a small energy gap of about 2meV at (0.5,0), which is 
an indication of single-ion anisotropy. Figures El^d) and 
[HJe) show the lower and upper limits of the gap in the 
upper band, and Fig. ^i) is a cut through the top of the 
dispersion at (1,0). 

We evaluated several different models by fitting the 
calculated magnon dispersion to that determined exper- 
imentally, shown in Fig. HJ^a) for the direction parallel to 
(/i, 0). To obtain the experimental dispersion we analysed 
the peak positions of approximately 50 cuts like those in 
Fig. [3] The spectrometer resolution was included in the 
analysis (see Supplemental Material [3^). 

The spin wave models are constructed from the Hcisen- 
berg Hamiltonian for the interaction between two spins 



at sites i and i-,T~i ~ 



• Sj . The spin Hamil- 



tonian for the Goodenough model has five exchange inter- 
actions Jij, four of which are in the Mn02 plane. These 
are shown in Fig.[lja), and are the nearest-neighbor (NN) 
FM interaction Jpi, the NN AFM interaction Ja, and the 
next-nearest-neighbor (NNN) FM interactions between 
two Mn"'+ sites Jf2 and between two Mn'^"'" sites Jps- 
There is also the NN AFM interaction between layers 
within a bilayer J± . Single-ion anisotropy is included by 
adding —D'^iiSf)'^ to H. Initially we will not distin- 
guish between Mn"^"*" and Mn^"*" and assume the spins on 
all Mn sites have the same value. The neutron scattering 
intensity was calculated by linear spin-wave theory. 

Spin wave dispersion curves computed from the best-fit 
Goodenough model are drawn in Fig.HJa), and Figs.[2][c) 
and[2l[d) show simulated intensity maps. Taking an aver- 
age Mn'^+ and Mn^"*" spin of = 7/4, we obtain fitted pa- 
rameters Jpi = 11.39±0.05meV, Ja = -1.50±0.02meV, 
Jf2 = 1.35 ± 0.07 meV, Jpg = -1.50 ± 0.05 meV, Ji_ = 
-0.88 ± 0.03 meV and D = 0.074 ±0.001 meV. As there 
are many spins in the magnetic unit cell, there are many 
branches in the dispersion diagram. To aid the interpre- 
tation of Fig. m^a) , we only plot modes whose calculated 
intensities are >10% of the maximum intensity. 

The striking feature of this fit is the dominance of 
the NN FM interaction Jpi, almost an order of mag- 
nitude larger than any other exchange. It is stronger 
than the FM exchange in metallic FM manganites, 
for which Jpi « 3.5-5.5 me V [2^, [l^, and similar to 
Jpi = 9.98 meV reported for the AFM CO phase of 
Lao.5Sri.5Mn04 [i^. Another notable result is that the 
NNN interactions between FM aligned spins (Jp2 and 
Jps) are almost identical in magnitude but opposite in 
sign, a feature imposed by the energy of the center of the 
gap and its size. This suggests some frustration in the 
ground state. The dispersion is unaltered if the values of 
Jp2 and Jp3 are exchanged. However, a comparison of 
the measured intensity with simulations in {h, k) planes 
at 60 and 35meV [see Figs. |li;b)-lli;d) and Figs. \MS)- 
m^h), respectively] unambiguously determines that Jp2 is 
FM and Jp3 is AFM. 

In the model just described, the gap in the upper band 
is caused by NNN interactions. Another possibility is 
that it could arise from alternating spin values due to 
charge disproportionation. Assuming that the number of 
d electrons is proportional to the total spin, with S ^ 2 
associated with Mn^"*" and 5 = 3/2 with Mn^+, we have 
found that charge states of Mn2-68(2)+ and Mn''-^^^^)-^ 
are required to reproduce the observed gap. In contrast, 
bond valence sums yield Mn^-^^(^)+ and Mn^-^^(^)+ 
a disproportionation only 25% of that required 
possibility, therefore, can be dismissed. 

A similar analysis was performed for the ZP model, 
except the spins in the Hamiltonian refer to the Zener 
polaron, rather than an individual Mn ion. We used 
an effective Heiscnberg Hamiltonian to capture the es- 
sential features of the ZP model [24|. It has only three 
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FIG. 4. (Color online), (a) Comparison of the measured magnon dispersion (symbols) parallel to the {h,0) direction with 
spin-wave calculations for the Goodenough model, (b)-(i) Measured and simulated intensity on {h, k) planes at the upper 
[(b)-(e)] and lower [(f )-(i)] limits of the gap. The intensity is averaged over 58-62 meV for the upper limit and 34-36 meV for 
the lower limit. Panels (b)-(e) show data, best-fit Goodenough model, best-fit Goodenough model with Jf2 and Jps exchanged, 
and best-fit dimer model, respectively. Panels (f)-(i) show the same but for the lower limit. Intensity units are mbsr~^ meV~^ 
per Mn. 



in-plane magnetic interactions [see Fig. IHb)]: the NN 
AFM interaction between two polarons (Ja), and two 
NN interactions reflecting the zig-zag structure ( Jfu and 
Jfd)- The out-of-plane exchange J± is also included. 
The magnetic ground state is CE-type (i.e. the ZP mo- 
ments are parallel along a zigzag) if Jpu > Jfb and 
\Ja\/Jf\j > Jfi)/{Jf\j — Jfb), and the ground state has 
the 90° structure shown in Fig. [IJb) when Jpu = Jfd *C 
I Ja|. These correspond to the $ = and $ = 7r/2 ground 
states of Ref. [l6| , respectively. In neither case can val- 
ues of the exchange parameters be chosen to to reproduce 
the observed 35-55 meV gap or dispersion maximum at 
h = 1. To illustrate this, Fig. [5] shows the spectrum 
of the best-fit Goodenough model alongside representa- 
tive simulations for the two ZP models. In the CE-type 
ZP order, the dispersion parallel to the zigzags has max- 
ima at h ~ 1/2,3/2,..., i.e., inconsistent with the data 
above the gap. In the 90° ZP phase, the quasi-lD dis- 
persion is perpendicular to the FM zigzags with maxima 
atfc = l/4,3/4,..., which because of twinning would also 
occur at these values of h, also inconsistent with the data. 
Furthermore, the weaker dispersion measured perpendic- 
ular to the zigzags confirms that Ja is not the dominant 
exchange parameter. 

A third model was tested, which is like the Goode- 
nough model but which allows for dimerization within 
the zig-zag FM chains. Mn ions are paired similar to 
the ZP model, but the model does not enforce the rigid 
coupling between the spin pairs. This dimer model is 
shown in Fig. [ijc) . The important interactions are the 
two NN FM interactions Jps and Jpw, the NN AFM in- 
teractions Jai and Ja2, and the intrabilayer interaction 
J±. As before, we assume all Mn spins are S = 7/4, and 
we set Jai = Ja2 = Ja- The fit of this model to the 



measured dispersion differs from the Goodenough model 
fit by ^ 2% and therefore is practically indistinguish- 
able. The fitted parameters arc Jps — 14.2 ± 0.08 meV, 
Jfw = 8.43 ± 0.06 meV, Ja = -1.52 ± 0.01 mcV, J± = 
-0.92 ± 0.03 meV, and D ^ 0.073 ± 0.001 meV. In the 
dimer model the gap arises from the alternating NN FM 
exchange without the need for NNN interactions. The 
dimer model is superficially appealing because of this 
simplicity and because the value of the intra-dimer ex- 
change Jps appears not unrealistic for double exchange 
(see above). In fact, the CE-type ZP dispersion cor- 
responds to the limiting case of the dimer model as 
Jps oo. However, a comparison of the measured 
intensity with simulations in {h, k) planes at the upper 
and lower limits of the gap [see Figs. Hl^b) and|4|e), and 
Figs.mjf) anddji)] rules out the dimer model. 

We emphasize that our study can conclusively distin- 
guish between different models only because it includes 
data above the 35-55 meV gap. This is demonstrated in 
Fig. [5l in which the spectrum for the dimer model in the 
strongly bound dimer limit Jpg 3> Jpw (equivalent to 
the CE-type ZP model) has been calculated with param- 
eters adjusted to match it to the lower spin-wave branch 
of the Goodenough model. Either model would be ac- 
ceptable if the data were restricted to energies below the 
gap, but only the Goodenough model is consistent with 
the branches observed above 40 meV. 

This study of the spin wave spectrum of PCSMO 
clearly demonstrates that the Zener polaron model is not 
a suitable description of the ground state in the AFM CO 
phase. The dimer model, which does not impose a rigid 
coupling between the two Mn ions in a Zener polaron, 
similarly fails to describe the data. What distinguishes 
our data from previously published data on AFM CO 
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FIG. 5. (Color online). Simulated spectra for the best-fit 
Goodenough model, the CE-type ZP model (which is equiv- 
alent to the strongly-bound dimer model), and the 90° ZP 
model. 



phase manganites [25|, |26| is that they cover the entire 
one-magnon spectrum, which was essential to rule out 
the ZP and dinier models conclusively. In contrast, the 
results are very well described with the original Goode- 
nough model, even without any significant charge order- 
ing. The entire magnetic spectrum can be explained by 
a single value of the NN FM interactions between Mn 
sites along the FM zigzags which is larger than that typ- 
ically found in metallic manganites, strongly suggesting 
its origin in double exchange. Significant NNN interac- 
tions along the zigzags further suggest the Cg electrons 
are not fully localised. 
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